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Coulomb correction to the screening angular parameter of the Moliere multiple scattering theory is 
found. Numerical calculations are presented in the range of nuclear charge 2 < Z < 82. Comparison 
with the Moliere result for the screening angle reveals up to 30% deviation from it for sufficiently 
heavy elements of the target material. 

PACS numbers: 

<N 

o 

^ ■ I. INTRODUCTION 

o 

^"^ ' The Coulomb correction (CC) is the difference between the exact in the parameter £ = Za//3 

result and the Born approximation. At intermediate energies, when kinetic energy of the scattered 
particles is approximately from 0.1 to 2.0 MeV, formulas for the Coulomb corrections are not available 
in analytical form [I]. The analytic formulas for the high-energy CC are known as the Bethe-Bloch 
'"^ , formulas for the ionization losses and those for the Bethe-Heitler cross section of bremsstrahlung 

A similar expression was found for the total cross section of the Coulomb interaction of compact 
| hadronic atoms with ordinary target atoms [91. Were also obtained CC to the cross sections of 

the pair production in nuclear collisions 0, Hull an interaction potential and the spectrum of 
bremsstrahlung 0, 0, • The specificity of the expressions presented in this work is that they define 
the CC to the Born screening angle Xa an d an exponential part of distribution function of the Moliere 
multiple scattering theory. 



(N 



The Moliere theory of multiple scattering [ll| is the most used tool for taken into account the 



multiple scattering effects in experimental data processing. The experiment DIRAC and many others 
[12-16] face the problem of excluding the multiple scattering effects in matter from obtained data. As 
the Moliere theory is currently used in the energy range roughly from 1 MeV to 200 GeV, the role of 
| the high-energy CC to the parameters of this theory becomes significant. 

Of especial importance is the Coulomb correction to the screening angular parameter, as this 
parameter also enters other important quantities in the Moliere theory. In his original paper, Moliere 
received an approximate semianalytical relation for the exact \ a an d first-order \a values of the 
r \ . screening angle: 

3: 



X a = xf^l + 3.34(Za//3) 2 



While the first term of this expression is defined quite accurately, the coefficient in the second term is 
found only numerically and approximately. 

In this work, we have obtained for \ a an exact in £ compact analytical result. In the second-order 
in £, it is given by: 



X a = Xa\/l + 2.l3(Za/pf 



We have also evaluated numerically, in the range 2 < Z < 82, the Coulomb corrections to the Born 
approximations of the screening angular parameter. Additionally, we have estimated the absolute and 
relative accuracies of the Moliere theory in determining these corrections. 

The paper is organized as follows. In Section 2, we consider the standard approach to the multiple 
scattering theory proposed by Moliere. In Section 3, we obtain the analytical and numerical results 
for the Coulomb corrections to the screening angular parameter. In Conclusion, we briefly summarize 
our results. 
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II. MOLIERE MULTIPLE SCATTERING THEORY 

Multiple scattering of a charged high-energy particle on the atoms of a target is the diffusion process 
in the angular plane of (6, <j>) = X- The angular phase volume is dO = 2nxdx, we suppose % = |%| <C 1. 
We define a(x)xdxd4> as the differential scattering cross section into the angular interval XiX + dx\ 
the Rutherford scattering cross section reads 

. . {2Ze 2 \ 2 1 JT e 2 Z 2 A 

<? R (X) = J — , U = , e 2 = Ana , 1 

\ rmr J x r 

where m and v are the mass of the scattered particle and its velocity at large distances from the 
scattering center, which assumed to be at rest; a is the fine structure constant and Z is the atomic 
charge number. 

Define now W(9, t)9d9 as the number of electrons scattered in the angular interval d9 after traveling 
through the target of thickness t. The normalization condition J W(9,t)d 2 9 = 1. The Boltzmann 
transport equation is 

' !U /:) = -n W(9,t) [ a( X )xdx + n f W{6 - X ,t)<r(x)d 2 x , (2) 



dt 



in which no is the number of the scattered centrum in 1 cm 3 , d 2 x = x^X^0/(2tt)- The first term in the 
right-hand side describes the decreasing in the number of electrons from the cone 9, and the second 
the increasing in the cone from the outside of the cone. 

Following the Moliere, we introduce the Bessel transformation of distribution 

oo 

g( V ,t) = J 9J (n6)W(9,t)d6 , (3) 
o 

oo 

W(9,t) = J vMv6)g(v,t)d V ■ (4) 



For g(r},i), using the folding theorem, we obtain 

oo 

dg(ii,t) 



dt 



= -n g(v,t) J o-(x)x^x[i - Mvx)} ■ 



(5) 



Its solution is 



.g(7?,t)-cxp{iV(7 7 )-iV }, (6) 
JVfa) = J a( X )xdxJo(vx) ■ (7) 

Inserting this expression in the Bessel transform of the distribution function W we have: 

oo ( oo *\ 

W{9, t) = ^ J vdvMvO) cxp j -n 1 2tt J a( X )xdx [1 - Mvx)] \ ■ (8) 
o I o J 

For the screening potential, the differential scattering cross section reads 

4Z 2 e 4 _ h_ 

{vp) 2 ( X 2 +Xo) 2 ' X ° P a 

and the total cross becomes 



<X) = ,„,^ 2 , 2 , „2x 2 XO = — , (9) 



f z \ j 47ra 2 (Ze 2 ) 2 
a = 2nJ <r(x)xdx = /fc ^ 9 1 ■ (10) 



{hv) 2 
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In the case of a thin target t < i s , where l s — l/(na), the distribution function may be written as 

W{9, t) = nta{6) = — a{6) . (11) 

Here, TV is the number of scattering centers, S is the area of the thin target surface (Na(8) -C S), 
and W(6, t) — Na(6)/S is the probability of single scattering. 

For large values of x, the- cross section cr(x) ~ Vx 4 decreases rapidly. It is a complicated function 
for x ~ Xo wit h 

Xn = - , A = — , a = 0.885 a Z" 1/3 . (12) 
a mv 

Here, a is the Bohr radius of the particle, and a is the Fermi radius of the atom. 

For the reasonable thickness, the width of the multiple scattering distribution is very large compared 
with x - 

Let us write 



q( X ) 2 m + l)zV 
' — ' *c = ^n t — 



n t(r(x)xdx = 2 XcXdx — 7~ , Xc = ^n t — ^ , e = 4?ra , (13) 



where z is the charge of the scattered particle. q(x) is the ratio of actual to Rutherford scattering 
cross sections. We replace Z 2 — > Z(Z + 1) keeping in mind the scattering on atomic electrons. 

The physical meaning of x c can be understood from the requirement that the probability of 
scattering on the angles exceeding x c is unity: 

oo oo 

f , s 4irn n te A Z(Z+ 1) f dx 
2n t J daix) = 2 ^ > J -* = 1 . (14) 

Typically, x c /Xo = 100- The quantity q(x) is equal to unity for large values of x > X c ari d tends 
to zero at x = 0. The main x values belong to the region x ~ Xo- It contains deviation from the 
Rutherford formulae due to the effects of screening of atomic electrons and the Coulomb corrections 
arising from multiphoton exchanges between the scattered particle and the atomic nuclei. 
We obtain in terms of x c 



- In gfat) = 2x1 I %q(x)[l-Jo(xv)} ■ (15) 
X 



.2 





To estimate the value of integral, we introduce (following [lj], some quantity k from the region 

(Xo:Xc): 

Xo < k < x c ■ (16) 
In the region x < k, we can use 1 — Jo(xv) — (xv) 2 /4: 

fc k 

%q(x)[l-Jo(xr))} = jr, 2 [ — ?(x) • (17) 
X d 4 7 x 



In the region fc < x, we can put cj'(x) = 1' 

/ — [i- Mxv)} = ]v 2 h(k v ) 

J X 4 



OO oo 

h{x) = Aj^[l-Mt)] = ^[l-J Q (x)}+2jP 1 (t), 

X X 

oo oo 

2 J ?sMt) = \ji{x) + J jJo(t), x = k v . (18) 
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Here, we used 



Ji(t) - -Ji(x)+ / -J((t) 



OC 

\jx{x) + /*(■*>(*)- 



Jl(t) 



, tj'^ — t,J() — J\ . 



Using J (t)dt/t = \n(2/x) - C E + 0(x 2 ) at x < 1 we obtain 



§[1 - Jo OCT)] - ^ 2 [l - C £ - ln(faj) + C)(i /.■„)- 



(19) 



with the Euler constant C E = 0.57721 . . . 

Considering the contribution of the region \ < k Moliere introduce the notion of the screening angle 

Xa • 



ln X a = I™ 

k— >oo 



L0 



— q(x) + ^ - ln/s 
X 2 



Hence, we obtain 



]ng(r),t)=N -N(n) = -{ Xc nf 



Introducing a new variable y = x c V we § e t 



iV - N( V ) = -y 2 



b - In ( -y' 



^"■| + 1 - 2C - Sl "(^ 



The Moliere transformed equation is 



(20) 



(21) 



(22) 



W{6,t)dd6 = XdX I ydyj {\y)exp <{ -^y 2 



-/. ■[ -In ( -y 2 



a = e/xc 



(23) 



This rather simple formula permits one us to develop an iteration procedure for W . Really, putting 
b = B — hxB and introducing the variables x = A 2 / B and u = y\f~B one can obtain the expansion of 
the distribution function in a power series in 1/B: 



W(6)6d6 = -4—9d9 
X Z C B 



with 



W w (x) + -W w {x) + X;W {2) {x) + . . 
B B z 



W {n \x) = — uduJoiuyG) 



i 2 In 



V(X^), 9 2 =x 2 c B 



The result of numerical integration (x) was obtained in papers by Moliere, Bethe and Scott lllll7| 
(see also [18|). In practice, the value of B as a solution of the transcendental equation b = B — ln£> 
is large enough B ss 5 to provide the convergence of the expansion series. 
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Let us investigate the expression ln g = N — Nq(i]): 
In g 



y_ 

4 



B - In B - In ■ 



y 2 B 



l 

B 



y 2 B ln y 2 B 



(24) 



Its minimal value is \ng n 



e corresponding to the value j/q/4 = e . Since e ss (Xc/Xa) i s °f 



the same order of magnitude as the number of collisions No, the accuracy of the final result of Moliere 
increases with the number of collisions as exp{— No/e}. As was shown by Bethe, 



XT: 



= AnNA — 



1-167x2 

6680 1 (Z + l)Z 1 / 3 z 2 
f3 2 A(l + 3.34 £ 2 ) 



(Z+l)Z 1 / 3 z 2 0.885 2 
f3 2 A 1.167 (1.13 + 3. 76£ 2 ) 



(25) 



Here, 6 « 8.8, A^ A = 6.02 x 10 23 cm 3 is the Avogadro number, A is the atomic weight, £ = Za/ [5 is 
the 'Born parameter', /3 = v/c, and the thickness t of the target is measured in units grammxcm 2 . 

Some comments. The quantity e b — (x c /x'a) 2 depends on the screening angle \' a = l-080% a . The 
screening angle incorporates the deviation of the potential from the Coulomb one from both the 
screening effect of atomic electrons and the Coulomb corrections. 

Estimating the scattering phases within the Thomas-Fermi model 

V(r)=± A(-) , a = Z~ 1/3 x 0.466 x 10 _8 cm , 

r \aJ 

3 

A(r') = 5>e-K 



ai = 0.1 , a2 = 0.55 , 03 = 0.35 , 
b x = 6 , b 2 = 1.2 , 63 = 0.3 , 



in the limit a — > 00 or the limit of small £, one obtains 



<?(x) 



Xo 



K h b 2 + (x/XoY 



(26) 



(27) 



In order to obtain a result valid for large £ and also for large angles x^ Moliere used the WKB 
method and rather a rough approximation in describing the Coulomb corrections. His result for the 
screening angle turns out to be only numerical and approximate. As a result it leads to the first term 
Xo x 1-13 of the expression for screening angle: 



X^=Xo(l-13 + 3.76C 2 ) 



(28) 



In the next section, we will use for these purposes the eikonal approximation and obtain an exact 
expression describing the Coulomb correction to the value of the Born-approximation screening angle 
Xl = \/TT3Xo- 



III. COULOMB CORRECTION TO THE SCREENING ANGULAR PARAMETER 

Remind now the relations for scattering amplitude in eikonal approximation: 

/(q) = ^~ Jd 2 b exp{-*qb}S(b) , 5(b) = exp{^(b)} - 1 , 

OO 

cj)(b) = - I dzU{b,z) , U(b,z) = Za/r , r = \/b 2 + z 2 , (29) 



where (j)(h) is the eikonal phase (see [19(, Appendix E). 



G 



It's convenient to define an interaction potential in the Landau-Pomeranchuk-Migdal effect theory 
(see Appendix A in Q): 

V(b) = nj (l-cxp{*qb})|/(q)| 2 A, \f(q)\ 2 d 2 q = da(q) . 

In the case of the screened Coulomb potential, we have the following expression for the eikonal phase: 



0(b) 



Za f 1 f r 1 Za 

/ dz — exp < — > = 2 Kq 

v J r la) v 



(30) 



where a is the Thomas-Fermi atom radius, and Ko(b/a) is the modified Bessel function. 

The equation for the potential V(b) can be written (after performing the angular integration) as 



Comparing this result with 



27m 



No - N(ti) 
n t 



[1 - J (qb))da(q) 



= / [1 - Mvx)}di( x ) , 



(31) 



(32) 



in which Nq — N(rj) = — In 17(77), we obtain the similarity when accept qb — rjx, q = pf], b = x/p> V — 
mv. 

So the problem of deviation of the potential V(h) from the Born one 



AV(h) 



,[y(b)] 



exp 



{#(b + x) - </>(x)]} - 1 + - [#b + x) - 4>{x)\ 



with A cc [V(b)] = V(h) - V B (b) is similar to our problem of deviation of the screening angle in the 
eikonal approximation from its Born value: 



ln^)] = A cc [lng( V )] = ^(xcvf&cc [in {x' a f 



(XcVf 



1 

2^ 



d A x 



(x + b) J 



2 \ ^ 



(33) 



with the Coulomb corrections A cc \lng(rj)] = \ng(r]) — \ng B (r)) and A cc [ln (xL)] = m (Xa) — m (xL) 
The two-dimensional integral calculated in [gj turns out to be 



/(0 = e 2 E 



£ = Za/P 



From 



we obtain 



or, equivalently, 



J n(n 2 + £ 2 ) 

Acc [In (Xa)]=/(0 
A cc [ln(^)] = Re[V'(l+^)] +C B , 



(34) 



(35) 



(36) 



with C E = —ip(l). Here, we use the smallness of the ratio x/a <C 1, the b ~ x <C a and apply 
the relevant asymptotic of the Bessel function Kq(z) — C — ln(z/2) + 0(z 2 ). The main reason of 
such derivation of relations (f33|) and ([36]) is the significantly different regions of contributions of the 
screening effects and Coulomb corrections. Really, the last ones play he main role in the region of 
small impact parameters, where the number of atom electrons is small and the screening effects are 
negligible. These results may also be obtained with by using the technique developed in [5|. 
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Finally, we get 

Xl = 1.13x§ 
which can be compared with the Moliere one 

Xl = 1.13x2 



1 + 2.131 



T 



1 + 3.341^ 



(37) 



(38) 



(2) (2) 

To calculate the second-order relative corrections to the first-order results, S cc and S M , which 
correspond to Eqs. (f3"T)) and (|38|) . respectively, and to investigate their Z-dependence, we first present 
these equations in the approximate form: 



(xL) « (xI) J 



1 + 1.204 



Za 

T 



(39) 



(Xa) ^ (Xa) 



1 + 1.670 




Then, (SSI and 0D]) become 



x(2) 
Una 



(X'a) 







=s 1.204 


f Za 


{X'a 


r 


\T 


aL 2) ( 




=a 1.670 


' ' Za 




r 






(40) 

(41) 
(42) 



(2) 

In order to compare our results and those of Moliere, we also define the absolute A CCM and relative 



S'ccm differences between the values of 8^ (x'a) an d ^cc(x'a)'- 



c(2) _ 



A 



(2) 



f(2) 



f(2) 



r(2) 



^(2) 
7(21 



(43) 



Table 1 presents the Z-dependence of the second-order corrections to the first-order results for 
j3 = 1. It shows that the values of the relative corrections (J^c f° r large-Z targets (Z ~ 80) does reach 
40%. Hence, it is also obvious that with the rise in the nuclear charge the absolute accuracy /X^cm 
of the Moliere theory in determining the relative CC to xL due *° the difference of the coefficients 
in Pdj) and (pt2"|) increases to approximately 16%, and the corresponding relative error 5 CCM does not 
depend on Z and is about 28%. During our analysis, we omit systematically the contribution of order 
a compared with ones of order 1. 

We can also calculate the exact in £ absolute correction A cc [ln (x'a)} = /(£) an d relative correction 
S cc(Xa) to trie Bo ™ screening angle (xq) B , 



&Cc(Xa) = } \b = ; ,\b = ex P [/(£)] - 1 = 



(Xa) (Xa 

as well as compare S cc with the Moliere result S M . 



(44) 



c(2) 



r(2) 



f(2) 



(2) 



(45) 



For this purpose, we must first calculate the values of the function /(£) = Rep0(l + + C E . The 
digamma series 



^(1 + = 1-^- 



£(-l)«[ C (n-l)]f 



(46) 
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where £ is the Riemann zeta function and |£| < 1, leads to the corresponding power series for Re ['0(1 + 
»0] =Re^(«)] and |£| < 2: 

1 OO 

Re[«] = 1 - C E - —— + [C(2n + f )]£ 2 ", (47) 

^ n=l 

and the function /(£) = £ 2 53^Li[ n ( n2 + £ 2 )] -1 can be represented at |£| < 2 as [20] 

1 OO 

/(0 = i- TT T2+E(- 1 ) n+1 [c( 2n+1 )]^ n ( 4g ) 

n— 1 

= 1 3_ + 0.2021 £ 2 - 0.0369 £ 4 + 0.0083 f - ... 

f + £ 2 

An equivalent way to estimate /(£) to four decimal figures is to present the sum from (|34[) in the 
following form 0|: 

OO 

E = (i+a" 1 + E(-^)" _1 [c( 2 «+ 1 )-i] ! ( 49 ) 



= (! + C 2 ) * + 0.20206 - 0.0369£ 2 + 0.0083£ 4 - 0.002£ 6 . 

Equation ([49]) is sufficient to evaluate this sum up to £ < 2/3 = 0.667. 

The calculation results for (|4"9"]l . function /(£) (|48p. the relative Coulomb correction (S cc (14"4"]| . 

_____ ______ 

their difference A CCM , and relative difference S CCM (|45ll with the Moliere correction S M (|4"2"]l at /? = 1 
are given in Table 2. 

It will be seen from Table 2 that for the light elements up to Z = 28, all the exact corrections coincide 
with the corresponding second-order corrections of Table 1. Beginning with Z — 42 the relative CC are 
lower than the above-mentioned Sac (Xo) > ^cc (Xa) i an d this discrepancy increases approximately to 
10% for the heavy elements with Z ~ 80. The ma gnitude of A cc [in (xL)] = f( Za ) is about 30 % for 
Z ~ 80. The size of the corresponding relative CC for these values of Z is approximately 40%. The 
absolute and relative differences with Moliere corrections increases to 20% and 34%, respectively, at 
Z ~ 80. 

Thus, in the case of scattering on large- Z targets such corrections to the Moliere result (|42[) as A CCM 
and S CCM become significant and should be taken into account in the description of experiments with 
nuclear targets. 



IV. CONCLUSION 

We have calculated the Coulomb correction A cc [in (xL)] an< i the relative Coulomb correction 
S C c{xL) to the screening angle x' a both analytically and numerically in the range 2 < Z < 82. We 
have found that these corrections are the order of 30% to 40% for Z ~ 80. Additionally, we evaluated 
the difference and relative difference between our results in determining the CC to x' a an d those of 
Moliere and found that they are about 20% and 30%, respectively, for heavy atoms of the target 
material. 
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Table 1. The Z-dependence of the second-order corrections defined by Eqs. 



ivl 


Zj 


Zj (1 


r(2) 
Occ 


r(2) 
Om 


A (2) 

U CCM 


r(2) 

"CCM 


Be 


4 


0.029 


0.001 


0.001 


0.000 


0.286 


Al 


13 


0.094 


0.011 


0.015 


0.004 


0.280 


Ti 


22 


0.160 


0.031 


0.043 


0.012 


0.280 


Ni 


28 


0.204 


0.050 


0.070 


0.020 


0.286 


Mo 


42 


0.307 


0.113 


0.157 


0.044 


0.280 


Sn 


50 


0.365 


0.160 


0.222 


0.062 


0.279 


Ta 


73 


0.533 


0.342 


0.474 


0.132 


0.278 


Pt 


78 


0.569 


0.390 


0.541 


0.150 


0.279 


Au 


79 


0.577 


0.400 


0.554 


0.154 


0.278 


Pb 


82 


0.598 


0.431 


0.598 


0.166 


0.279 



Table 2. The Z-dependence of the exact with respect to £ corrections 

(133, 61, and (05]) 



M 


Z 


Za 


S 


f(Za) 


Sec 


^CCM 


&CCM 


Be 


4 


0.029 


1.201 


0.001 


0.001 


0.000 


0.286 


Al 


13 


0.094 


1.193 


0.011 


0.011 


0.004 


0.280 


Ti 


22 


0.160 


1.176 


0.031 


0.031 


0.012 


0.280 


Ni 


28 


0.204 


1.160 


0.049 


0.050 


0.020 


0.287 


Mo 


42 


0.307 


1.113 


0.105 


0.110 


0.047 


0.297 


Sn 


50 


0.365 


1.080 


0.144 


0.154 


0.068 


0.306 


Ta 


73 


0.533 


0.971 


0.276 


0.318 


0.157 


0.330 


Pt 


78 


0.569 


0.947 


0.307 


0.359 


0.182 


0.336 


An 


79 


0.577 


0.941 


0.313 


0.367 


0.187 


0.337 


Pb 


82 


0.598 


0.926 


0.332 


0.393 


0.205 


0.342 



